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Abstract 

All systems in thermal equilibrium exhibit a spatially variable energy landscape due 
to thermal fluctuations. Thus at any instant there is naturally a thermodynamically 
driven localization of energy in parts of the system relative to other parts of the system. 
The specific characteristics of the spatial landscape such as, for example, the energy 
variance, depend on the thermodynamic properties of the system and vary from one 
system to another. The temporal persistence of a given energy landscape, that is. 



'Permanent address: Departament de Qui'mica-Fisica, Universitat de Barcelona, Avda. Diagonal 647, 
08028 Barcelona, Spain 

^Present address: Max-Planck Institut fiir Festkorperforschung, Heisenbergstr. 1, 70569 Stuttgart, 
Germany 

■'■Permanent address: Institute de Astronomi'a, Apdo. Postal 70-264, Ciudad Universitaria, Mexico D. F. 
04510, Mexico 



1 



the way in which energy fluctuations (high or low) decay toward the thermal mean, 
depends on the dynamical features of the system. We discuss the spatial and temporal 
characteristics of spontaneous energy localization in ID anharmonic chains in thermal 
equilibrium. 



1 Introduction 

The pioneering work of Fermi, Pasta and UlamS demonstrated that a periodic lattice of 
coupled nonlinear oscillators is not ergodic, and that energy in such a lattice may never 
be distributed uniformly. A great deal of work has followed that classic paper trying to 
understand how energy is distributed in discrete nonlinear systems. M Specifically, the 
possibility of spontaneous energy localization in perfect anharmonic lattices has been a 
subject of intense interest. The existence of solitons and more generally of breathers 
and other energy- focusing mechanisms, and the stationarity or periodic recurrence or even 
slow relaxation of such spatially localized excitations, are viewed as nonlinear phenomena 
with important consequences in many physical systems0llllil 

The interest in the distribution and motion of energy in perfect arrays arises in part 
because localized energy in these systems may be mobile, in contrast with systems where 
energy localization occurs through disorder. The interest also arises because such arrays may 
themselves serve as models for a heat bath for other systems connected to them.0 Albeit 
in different contexts, "perfect" arrays serving as energy storage and transfer assemblies for 
chemical or photochemical processes are not uncommon. 

The study of anharmonic chains and of higher-dimensional discrete arrays has been 
less than systematic, certainly an inevitable consequence of the breadth and mathematical 
difficulty of the subject. Some studies (including the work of Fermi, Pasta, and Ulam) deal 
with micro canonical arrays. Here one observes the way in which a given constant amount 
of energy distributes itself among the elements of the array. The notion of "temperature" 
usually does not enter in these discussions, although such an association could be made if 
the energy is randomly distributed. Other studies of anharmonic chains (far more limited 
in number) deal with systems subject to external noise and other external forces. The 
questions of interest here involve the ways in which noise can enhance (as in noise-enhanced 
signal propagationlli) or even totally modify (as in noise-induced phase transitionJl^) the 
properties of the nonlinear array. Even more limited has been the study of systems that 
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are in thermal contact with one or more external heat baths maintained at a constant 
temperature.0lli Here the questions usually revolve around the robustness against thermal 
fluctuations of stationary or quasi-stationary solutions of the microcanonical problem. In 
both microcanonical and canonical systems, some work concentrates on stationary states 
or long-time behavior or equilibrium properties of the array, while other work deals with 
transport properties or with the approach to equilibrium. Furthermore, there is variation 
in the portion of the potential where the nonlinearity resides. Thus, in some cases the 
elements of the array are themselves nonlinear while in others it is the coupling between 
elements that is nonlinear (and, on occasion, both are nonlinear). 

Within this broad setting, our interest in this paper focuses on one-dimensional arrays 
of classical oscillators in thermal equilibrium^ An understanding of thermal equilibrium 
properties and the effects of nonlinearities on these properties is a prerequisite to the per- 
haps more interesting analysis of the nonequilibrium behavior of anharmonic lattices in 
the presence of thermal fluctuations and the approach to equilibrium in such systems. In 
particular, here we deal with the case of "diagonal anharmonicity," that is, the nonlinear- 
ity in our model is inherent within each oscillator in the array (representing, for example, 
intramolecular interactions) , while the connections between oscillators (representing, for ex- 
ample, intermolecular interactions) are ordinary linear springs. The anharmonicity may be 
soft or hard. We explore the conditions that lead to spontaneous energy localization in one 
or a few of the oscillators in the array, and the time it takes for a given energy landscape 
to change to a different landscape. One could undertake a parallel study in systems with 
anharmonic interactions between oscillators ("off-diagonal" anharmonicity). We address 
such systems in subsequent work.0 

The energy landscape is determined by the local potential of each oscillator, and by 
the channels of energy exchange in and out of each of the oscillators. The couplings be- 
tween oscillators provide one such exchange channel, and the coupling of the array with 
the heat bath provides the other. We shall see that different arrays (soft, hard) behave 
very differently in response to these channels. We broadly anticipate our conclusions by 
revealing that 1) persistent energy localization occurs in arrays of weakly coupled soft os- 
cillators even when strongly coupled to a heat bath (while such localization is absent in the 
hard chain); 2) persistent localization occurs in strongly coupled hard arrays provided they 
are weakly coupled to a heat bath (while such localization is absent in the soft chain); 3) 
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quasi-dispersionless mobility of localized energy requires off-diagonal anharmonicity. 

These remarks point to the fact that our analysis of anharmonic chains in thermal 
equilibrium could start from two "opposite" viewpoints. On the one hand, we might start 
by analyzing uncoupled oscillators in thermal equilibrium and then proceed to investigate 
what happens if we couple these oscillators to one another. This approach focuses on the 
entropic localization mechanismlll and the way in which the coupling between the oscillators 
eventually degrades it. On the other hand, we might start with a coupled isolated chain, 
focus on energetic localization mechanisms in such a cham0 and then proceed to investigate 
the ways in which thermal fluctuations and dissipation affect such local structures. Since 
we are explicitly interested in localization in the presence of thermal fluctuations, and since 
entropic effects have received far less attention than energetic ones, we choose to follow the 
former approach. 

No matter the sequence of our queries, since here our interest lies mainly in understand- 
ing energy localization in a nonlinear discrete array in thermal equilibrium and the way in 
which thermal effects depend on system parameters, we pose our questions as follows: 

• How is the energy distributed in an equilibrium nonlinear chain at any given instant 
of time, and how does this distribution depend on the anharmonicity? In other words, 
can one talk about spontaneous energy localization in thermal equilibrium, and, if so, 
what are the mechanisms that lead to it? 

• How do local energy fluctuations in such an equilibrium array relax in a given oscil- 
lator? Are there circumstances in the equilibrium system wherein a given oscillator 
remains at a high level of excitation for a long time? 

• Can local high-energy fluctuations move in some nondispersive fashion along the ar- 
ray? In other words, can an array in thermal equilibrium transmit long-lived high- 
energy fluctuations (if indeed they exist) from one region of the array to another 
without too much energy loss to dispersion? 

The answers to these questions have not been found analytically, and are for that reason 
most clearly presented in comparative fashion. Starting with an ensemble of uncoupled 
oscillators at thermal equilibrium, one knows exactly the behavior of a single harmonic 
oscillator and can say a great deal about the behavior of a single anharmonic oscillator 
from general thermodynamic considerations. Thus, for instance, the mean energy of a single 
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harmonic oscillator in thermal equilibrium at temperature T is E = ksT {ks =Boltzmann's 
constant). This energy is on average divided equally between kinetic and potential (a 
partition that enters importantly in questions concerning landscape persistence). A simple 
virial analysis immediately shows that a soft anharmonic oscillator in thermal equilibrium 
has energy greater than ksT while a hard anharmonic oscillator has energy smaller than 
ksT. Both share the property of the harmonic oscillator that the average kinetic energy 
is k^T jl^ but their average potential energies differ. One also knows exactly the energy 
fluctuations in a harmonic oscillator: the energy variance o"^ is equal to k'^T'^, and the 
ratio of (T to is therefore independent of temperature. The energy fluctuations are easily 
determined to be greater in a soft oscillator and smaller in a hard oscillator. From these 
facts one can arrive at rather definitive qualitative conclusions regarding the distribution and 
persistence of energy in ensembles of single oscillators and the effects of the anharmonicities 
on these features .111 

The situation becomes more complicated when such oscillators are connected to one 
another. Not only can the oscillators now exchange energy with the heat bath, but there 
are also coupling channels whereby oscillators can exchange energy with one another. The 
interplay of these various energy exchange channels and the effects of anharmonicity on this 
interplay are some of the issues to be addressed in this work. 

This paper is organized as follows. In Section ^ we introduce our model and notation. 
We fix some of the parameter values and briefly discuss the numerical methods used in our 
simulations. Here we introduce the hard, harmonic, and soft local potentials to be compared. 
In Section ^ we review and illustrate previous results for uncoupled oscillators in thermal 
equilibrium so as to establish the background for the coupled systems. The phenomenon of 
"entropic localization," whereby ensembles of single thermalized soft oscillators localize and 
retain energy more effectively than harmonic or hard oscillators, is recalled. In Section Q 
we explore the consequences of coupling our oscillators. In Section ^ we briefly address the 
mobility of energy fluctuations in our systems. Finally, Section ^ summarizes our flndings 
and anticipates further studies. 
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2 The Model and Numerical Methods 



Our system is a one-dimensional chain of N identical unit-mass oscillators labeled i = 
1,2, ••• ,N with harmonic nearest-neighbor interactions and on-site potentials V{xi) that 
may be hard, harmonic or soft. Here Xi is the displacement of oscillator i from its equilibrium 
position, with associated momentum pi. We assume periodic boundary conditions. The 
Hamiltonian of the system is 

i=i ^ ^ 

where k is the intermolecular force constant. Figure || is a schematic of the model. 



Figure 1: Illustration of the ID chain considered in this work. Each oscillator in the chain 
experiences an on-site potential and is harmonically bound to its nearest neighbors. 

To represent the thermalization of our chain the model is further expanded to include the 
Langevin prescription for coupling a system to a heat bath at temperature T via fluctuating 
and dissipative terms. The stochastic equations of motion for the chain are then given by 
the Langevin equations 

Xi = -k(2xi - Xi+i - Xi-i) - -yxi - ^^^^'"^ + rjiit) (2) 

axi 

where a dot represents a derivative with respect to time. The r/j(i) are mutually uncorrelated 
zero-centered Gaussian (5-correlated fluctuations that satisfy the fluctuation-dissipation re- 
lation: 

{rii{t)) = 0, imitHit')) = 2jkBT6ij6{t - 1') . (3) 

Since we are interested in assessing the effects of anharmonicities on energy localization. 
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Figure 2: Left panel: the on-site potentials defined in Eq.(^). Right panel: the associated 
forces. Solid lines: harmonic potential, Vo(x). Dotted lines: soft anharmonic potential, 
Vs{x). Dashed lines: hard anharmonic potential, Vh{x). 

we start by specifying the on-site potentials to be used in our analysis: 




Vs{x) = \x\ - ln{l + \x\) (4) 
Vh(x) = -x^ + -x^. 

The subscript stands for the harmonic case, s for the soft and h for the hard. At small 
amplitudes the three potentials are harmonic with a unit natural frequency. Fig. ^ shows 
the potentials and associated forces. 

We end this section with a brief description of the numerical methods used in our 
simulations throughout this paper. The numerical integration of the stochastic equations for 
all our simulations is performed using the second order Heun's method (which is equivalent 
to a second order Runge Kutta integration) Ellil We use a time step At = 0.005. The 
number of oscillators in our simulations ranges between 100 and 1000 and is indicated in 
each figure as appropriate. In each simulation the system is initially allowed to relax for 
enough iterations to insure thermal equilibrium, after which we take our "measurements." 
In all of our subsequent energy landscape representations we have used the same sequence 
of random numbers to generate the thermal fluctuations. 
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3 Properties of Uncoupled Oscillators: Entropic Localization 



In order to understand the equilibrium properties of a chain of oscillators it is useful to first 
review the behavior of single (uncoupled) oscillators described by the potentials in Eq. 

Suppose first that our oscillator is isolated. The salient features of anharmonic oscillators 
are that 1) they oscillate with different frequencies at different energies, and 2) the density 
of states changes with changing energy. In particular, hard potentials are associated with 
increasing frequencies of oscillation and sparser densities of states with increasing amplitude 
(energy); on the other hand, in soft potentials the oscillation frequency decreases and the 
density of states increases with increasing amplitude. 
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Figure 3: Oscillation characteristics of single isolated oscillators. Left panel: frequency 
as a function of the oscillator energy for the potentials in Eq. (^. Right panel: oscillation 
periods for single oscillators. Solid lines: harmonic oscillator. Dotted lines: soft anharmonic 
oscillator. Dashed lines: hard anharmonic oscillator. 

To get a sense, useful for later analysis, of these and associated oscillator characteristics, 
we present several figures that show various distinct features of our three types of oscillators. 
Figure ^ shows the frequencies uj{E) of isolated single oscillators as a function of increasing 
energy E (which in turn corresponds to increasing amplitude). This frequency is evaluated 



directly by solving the equation of motion dx/dt = ±y/2[E — V{x)] over one period of 
oscillation at energy E: 

/ , . (5) 
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The amplitude of oscillation Xmax at a given energy can be found by solving for the positive 
root of V{x) = E. The harmonic oscillator has a single frequency at unity. The soft 
and hard oscillators oscillate at unit frequency at low amplitudes (energies) because we 
have chosen all the oscillators to coincide there, but with increasing amplitude the hard 
oscillator frequencies increase and those of the soft oscillator decrease. In Fig. ^ we also 
show the period of oscillations t{E) = 2'k/i^{E). The period increases with increasing 
energy for the soft oscillator, remains constant for the harmonic oscillator, and decreases 
with energy for a hard oscillator. This behavior will figure prominently in our subsequent 
analysis of energy localization. 

Next we consider these same single oscillators, but now each connected to a heat bath 
at temperature T via Langevin terms. The left panel of Fig. ^ shows the normalized energy 
distribution P{E) vs E for the three cases. This distribution is given by 

^^^^ " dEe-E/''BT^(^E^ (6) 

where the density of states is just the period of oscillations. The figure supports our in- 
troductory comments: firstly, that the average energy of the soft oscillators is greater than 
that of the harmonic oscillators, whose average energy is in turn greater than that of the 
hard oscillators; secondly, that the energy fluctuations are smallest in the hard oscillator 
and largest in the soft oscillator. Thus in equilibrium we find at any instant that there is a 
greater variability of energy in an ensemble of single soft oscillators than in one of harmonic 
or hard oscillators. The right panel of Fig. ^ shows the average period of oscillation T{kBT) 
for a thermalized distribution: 

POD 

T{kBT) = / dE t{E)P{E) . (7) 

JO 

Consonant with the energy dependence of t{E), the average period of the soft oscillator 
increases with temperature, that of the harmonic oscillator is independent of temperature, 
and that of the hard oscillator decreases with temperature. 

The features just discussed are also visible in the energy landscape rendition shown in 
Fig. ^. Along the horizontal direction in each panel lies an ensemble of 100 independent 
thermalized oscillators and the vertical upward progression shows how these oscillators 
evolve with time in the equilibrium system. Here and in all our energy landscape figures 
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Figure 4: Left panel: energy distribution in single thermalized oscillators for the three 
potentials at ksT = 0.5. Right panel: average oscillation period for the three oscillators 
as a function of temperature. Solid lines: harmonic potential, Vo{x). Dotted lines: soft 
anharmonic potential, Vs{x). Dashed lines: hard anharmonic potential, Vh{x). 

the y axis covers 120 time units, the same units shown on time axes throughout the paper. 
Each oscillator is connected to a heat bath. The grey scale represents the energy - an 
oscillator of higher energy is darker in this portrayal. 

The first thing to note is that along any horizontal line (i.e. at any given time) the 
soft landscape is darker and grainier than the harmonic, and the lightest and least grainy 
is the hard oscillator landscape. This reflects the fact that the soft oscillators have the 
highest energies and the greatest energy fluctuations. This observation provides a basis 
to be used in answer to the first question posed in the introduction. In an ensemble of 
independent oscillators in thermal equilibrium there is of course a greater energy in some 
oscillators than in others simply because there are energy fluctuations in a system in thermal 
equilibrium. These fluctuations are greater in soft anharmonic oscillators than in harmonic 
or hard anharmonic oscillators. 
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Figure 5: Energy (in grey scales) for ensembles of 100 thermalized independent oscillators 

as a function of time. The oscillators arc lined up (but not connected) along the x-axis 
and time advances along the y-axis. The temperature is ksT = 0.5 and the dissipation 
parameter is 7 = 1. Top panel: soft oscillators; middle panel: harmonic oscillators; lower 
panel: hard oscillators. 
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Figure 6: Energy correlation function vs time for independent oscillators with A:bT = 0.5 
and 7 = 1. Note that the energy changes most slowly in the soft potential ensemble. 
Solid line: harmonic potential. Dotted line: soft anharmonic potential. Dashed line: hard 
anharmonic potential. 

The second noteworthy feature of the landscape illustrates the answer to the second 
question posed in the introduction, namely, how long it takes in an equilibrium ensemble for 
the fluctuations to relax and the energy landscape to change. The trend for our independent 
oscillators is clear: the soft oscillators retain a given energy for a longer time than do the 
harmonic, which in turn hold on to a given energy longer than do the hard oscillators. This 
is particularly evident for those oscillators that acquire a high energy through a fluctuation: 
in the soft oscillator landscape the dark streaks are clearly visible. The reason for this 
behavior becomes clear if we write the equation of motion for the energy E = + V{x) 
for each oscillator. Setting p = x and using Eq. (||) one finds that for any type of oscillator 



Thus, the energy exchange with the surroundings involves only the momentum variable (i.e., 
the kinetic energy) . Consider an oscillator that has acquired a given high-energy fluctuation 
E, and consider how this energy is distributed between the oscillator displacement and 
momentum. In a harmonic oscillator the energy during one cycle of oscillation is equally 
partitioned between kinetic and potential. In a soft oscillator, however, the energy spends 



E = — +pri{t). 



(8) 
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relatively more time in potential than in kinetic form (and the opposite is true for the hard 
oscillator). Thus, during the major portion of the cycle the momentum of a soft oscillator is 
relatively low (while its displacement is large); the energy in the soft oscillator can therefore 
not enter from and leave to the thermal surroundings as easily as in the other oscillators. 
The energy relaxation process is therefore slower, and a soft anharmonic oscillator retains 
a high energy it might have gained via a fluctuation for a longer time.0 

The energy relaxation process visible in Fig. H is shown more quantitatively in Fig. |^. 
Here we have plotted the normalized energy correlation function 



The inner brackets indicate an average over time t (200,000 iterations) and the outer brackets 
an average over an ensemble of 1000 oscillators. The correlation function is normalized so 
that all energies, high and low, contribute "equally." It is thus a measure of the full exchange 
of energy with the heat bath, both through the dissipative term and also via the fluctuations. 



note that on average the energy of an oscillator changes on the time scale of half a period 
of oscillation, i.e. on the time scale it takes the oscillator to move from one side of the 
potential well to the other. 

We have thus summarized and illustrated our earlier findings ,0 namely, that in an 
array of independent oscillators in thermal equilibrium at a given temperature there are 
larger energy fluctuations and longer retention of energy the softer the oscillators. This 
is an entropy-driven localization, arising from the fact that the density of states in soft 
oscillators increases with increasing energy. It minimizes the free energy because it is 
entropically favorable for oscillators to populate phase space regions where the density of 
states is higher, which in an ensemble of soft oscillators leads to a greater spatial variability 
than in harmonic or hard oscillators. The temporal persistence of this greater variability is 
a consequence of the fact that coupling to a heat bath occurs only via the kinetic energy. In 
the soft ensemble the energy is in potential form a greater fraction of time than in kinetic 
form, which is not the case for the other ensembles. 




(9) 



We note that the trend in Fig. |6| (slower decay as the oscillators soften) is consistent with 
the corresponding slowing trend for each temperature in the right panel of Fig. ^. Also 
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Figure 7: Energy landscapes for thermalized independent soft oscillators as a function of 
time for different temperatures. The dissipation parameter is 7 = 1. Temperatures from 
top to bottom: fesT = 0.1, 0.5, 1.0 and 2.0. 
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We gave this scenario the name stochastic localization in our earher work,t3 but wiU 
refer to it as entropic localization, a term that more accurately reflects its physical causes. 
It is important to stress that entropic localization in soft oscillators is robust in the sense 
that it becomes more pronounced as temperature increases provided the potential continues 
to soften, and that it is achieved regardless of the initial condition of the system. 

The remaining parameters that can be varied at this point are the dissipation parameter 
and the temperature. A change in the dissipation parameter does not affect Fig. ^ since this 
is an equilibrium distribution. In Fig. [s] a higher dissipation parameter would cause a more 
rapid decay of energy fluctuations (and, correspondingly, a lower dissipation parameter 
allows an energy fluctuation to survive for a longer time). Thus, although high dissipation 
does not interfere with the appearance of greater energy fluctuations in the soft oscillators, 
it works against the temporal retention of excess energy by any one oscillator. The energy 
correlation function decays more slowly for the soft oscillator for any dissipation, and this 
decay is more rapid (for all the oscillators) as the dissipation increases. In any case, for a 
given dissipation parameter the softer potential retains energy for a longer time. 

The temperature affects the quantitative outcome of Figs. |^ and ^. In Fig. ^ the dis- 
tributions broaden with increasing temperature, but the differences between the different 
oscillators remain and, in particular, the fact that the distribution for the soft oscillator is 
the broadest continues to be true. In Fig. |5| higher temperatures produce relatively greater 
graininess in the soft oscillator figure than in the other two. This is clearly observed in 
the sequence of Fig. 0, which shows the evolution of ensembles of soft oscillators for differ- 
ent temperatures. A temperature increase leads to stronger entropic localization and this 
effect also appears in the energy correlation functions, as shown in Fig. ^. This behavior 
is contrasted with that of harmonic and hard anharmonic oscillators, whose energy land- 
scapes and energy correlation functions show essentially no temperature dependence in this 
range. The energy fiuctuations in these latter cases dissipate very quickly. Note that the 
temperature dependence of the correlation times implicit in Fig. ^ is consistent with the 
temperature dependence of an average period of oscillation of a soft oscillator as shown in 
the right panel in Fig. |^: with increasing temperature the correlation time continues to be 
approximately half a period. 

With this background, we are now ready to consider the behavior of chains of oscilla- 
tors, where everything that we have found so far has to be reconsidered in the face of the 
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additional forces now present through the oscillator-oscillator coupling. 

4 Coupled Oscillators 

In this section we explore the consequences of coupling the oscillators discussed in the 
previous section with harmonic springs. In this exploration we attempt to bring some order 
to seemingly contradictory reports that the coupled oscillators must be hard in order for such 
an array to localize energy effectively, or that the coupled oscillators must be soft in order 
to accomplish such localization. To anticipate our results: we will show that both claims 
are correct, but each in a different parameter regime and for different physical reasons. The 
variable parameters in this discussion are the temperature ksT, the dissipation parameter 
7, and the coupling strength k. 
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Figure 8: Energy correlation function vs time for independent soft oscillators with 7 = 1 
and different temperatures (the same as in Fig. ^. 

In order to determine the conditions that may lead to energy localization in a thermalized 
chain of oscillators it is useful to investigate the ways in which energy may escape from a 
given oscillator. It is apparent from the Langevin equation (^) that there are now two 
channels of escape. As in the last section, one is the friction term that dissipates the energy 
to the bath. The other is the coupling term that transfers energy to the nearest neighbors. 
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The difference between these two mechanisms is that the dissipation is determined entirely 
by the kinetic energy of the oscillator. Energy transfer along the chain, on the other hand, 
while still dependent on the kinetic energy, is primarily determined by the extension or 
contraction of the springs connecting neighboring oscillators, that is, by the potential energy 
through the relative oscillator displacements. To make these statements more quantitative, 
it is useful to generalize the concept of a local energy by defining a local function whose sum 
over sites is the total energy of the chain. To include the contribution from the nearest- 
neighbor restoring forces one writes 



and the total energy of the system is then E = ^^Ei- The rate of change of the local 
energy is easily found to be 



Note that although this expression does not explicitly involve the potential, the rate of local 
energy loss of course does depend on the potential through the displacements and momenta. 

The dynamics of the local energy will thus depend on the interplay of the thermal 
(fluctuations), dissipative, and intrachain forces. In order to highlight the main comparisons 
and contrasts, we frequently will juxtapose the behavior of chains for which one or the other 
of the energy exchange channels is clearly the dominant one, and in each case assess the 
effects of temperature changes. 

The effect of interoscillator coupling on entropic localization is illustrated in Fig. |9[ In 
this figure we show the system of soft oscillators that were uncoupled in Fig. |^ (specifically, 
the case with ksT = 0.5 and 7 = 1), but now providing successively larger values for the 
coupling constant k. Entropic localization is still apparent for small values of k, but as 
coupling increases there is clear degradation of entropic localization. This is to be expected 
since energy exchange is sensitive to large oscillator amplitude differences in soft oscillators. 

The associated energy correlation functions 
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Figure 9: Energy landscapes for thermalized soft oscillators as a function of time. The 
dissipation parameter is 7 = 1 and the temperature ksT = 0.5. Prom top to bottom the 
coupling constants are k = 0, 0.05, 0.5 and 1.0. 
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for the cases in Fig. |^ are shown in Fig. These curves confirm the degradation of entropic 
localization with increasing k. 

We thus turn to chains of coupled oscillators with low dissipation (7 = 0.05) and focus, 
in particular, on strongly coupled chains (if both k and 7 are small we know pretty much 
what happens from the analysis in the previous section). In Fig. |ll| we have drawn the 
energy landscape for the soft (top panel), harmonic (middle panel) and hard oscillators 
(lower panel) providing ksT = 0.5, 7 = 0.05, and k = 1.0. From this figure it is clearly 
evident that now the localization of energy at a given site is greater in the hard case than 
in the harmonic case, and this in turn, is greater than in the soft case. The confirming 



local energy correlation functions for these cases are plotted in Fig. |12|. Clearly, for a given 
temperature the hard array retains energy at a given location for a longer time than do the 
other two arrays. 



1.0 F 




Figure 10: Local energy correlation function vs time for chains of soft oscillators with 7 = 1, 
ksT = 0.5 and different values for the coupling constant (same as in Fig. ^). 

In the I0W-7, large-A: regime the effective energy exchange channel is sensitive to the 
oscillator amplitude rather than to its kinetic energy, so we expect entropic localization in 
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Figure 11: Energy landscapes for thermalized strongly coupled oscillators as a function of 
time. The dissipation parameter is 7 = 0.05, the temperature ksT = 0.5, and the coupling 
constant k = 1.0. Top panel: soft oscillators; middle panel: harmonic oscillators; lower 
panel: hard oscillators. 
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the soft array to be degraded since soft oscillators have large ampHtudes. Furthermore, 
as the harmonic couphng increases it eventually overwhelms the local soft potential and 
the soft chain becomes an essentially harmonic chain at sufficiently large k. On the other 
hand, hard oscillators exchange little energy via the coupling channel since they do not 
reach large amplitudes. This, and the fact that dissipation to the bath via kinetic energy 
(the other energy exchange channel) has been minimized (low 7), leads to persistent energy 
localization in the hard array. This is an energetic localization mechanism. The frequency 
mismatch between an energetic hard oscillator and its less energetic neighbors, and the 
dearth of density of states at high energies, further contribute to this persistence. 




Figure 12: Energy correlation function vs time for coupled oscillators with 7 = 0.05, ksT = 
0.5, and k = 1.0. Solid line: harmonic potential. Dotted line: soft anharmonic potential. 
Dashed line: hard anharmonic potential. 

The energetic localization mechanism in strongly coupled hard oscillators is robust 
against temperature increases. Indeed, according to our explanation, the localization should 
become more pronounced and persistent as temperature increases provided the dissipation 



is sufficiently weak. In Fig. 13 we have drawn the energy landscapes for a strongly coupled 
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Figure 13: Energy landseapes for thermalized coupled hard oscillators as a function of time. 
We take 7 = 0.05 and k = 1.0. Temperatures from top to bottom: ksT = 0.1, 0.5, 1.0 and 
2.0. 
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(k = 1.0) array of hard oscihators, weakly coupled to the bath (7 = 0.05) at different 
temperatures. The figure qualitatively confirms these expectations. The corresponding 
energy correlation functions are plotted in Fig. |14|: C(t) for the hard chain does decay 
more slowly with increasing temperature. Thus, localization in this strongly coupled system 
of hard oscillators becomes more effective with increasing temperature and is not entirely 
fragile against dissipative forces. On the other hand, the soft and harmonic correlation 
functions (not shown here) are essentially independent of temperature. Note that the trend 



in Figs. |12| and 14 is "opposite" to that of the uncoupled oscillators in the right hand panel 
of Fig. ^. In the strongly coupled chain harmonic and soft oscillators in fact lose their energy 
rather quickly on the time scale of one oscillation period of an isolated oscillator, but the 
hard oscillators retain energy correlations for longer than a period, indeed for many periods 
at the highest temperatures shown. With increasing temperature the hard oscillators retain 
energy more effectively even while the average oscillation period decreases. In fact, the 
decay of the correlation functions appears to involve two time scales, one of the order of an 
oscillation period and another much longer one that grows with temperature. 

The temporal irregularities (oscillations) visible in Figs. |l^ and 14 are reproducibly there 
at all temperatures; we do not know their source. 

5 Mobility of Localized Energy Fluctuations 



The first two energy landscapes in Fig. 11 show what might appear as fairly dispersionless 
energy transport. Narrow high-energy pulses move visibly along the chain before disappear- 
ing, while others appear (via thermal fluctuations) to repeat the process elsewhere along the 
chain. However, this can not be claimed to represent nonlinear behavior since the middle 



panel in Fig. 11 in fact represents completely harmonic system! This serves as a cautionary 
note about the overinterpretation of such results. 

We noted earlier that with increasing k the soft chain eventually becomes essentially 
harmonic because the intermolecular harmonic interactions overwhelm the local soft poten- 
tial (the only way to prevent this is by considering soft interoscillator interactions, which 
we defer to another pape r§). The upper panel in Fig. O exhibits mostly this essentially 
harmonic behavior - it is quite similar to the middle panel ~ but not entirely so. The 
soft oscillator chain clearly shows higher-energy regions than the harmonic (darker patches. 
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Figure 14: Energy correlation function vs time for strongly coupled hard oscillators with 
7 = 0.05, k = I and different temperatures (same as in Fig. [l3|) . 

a not fully degraded remnant of entropic localization) that move more rapidly (steeper 
streaks) over longer distances (longer streaks) than in the harmonic chain. Therefore, the 
soft anharmonicity is clearly still playing some role, albeit a diminishing one with increasing 
coupling. To provide some quantification, we introduce the dynamical energy correlation 
function 



, _ / {E.im+jjt + r)) - {E,it)){E,+jit + r)) \ 

^) - ^ (Em) - {Em? ^ • ^ 



This correlation function plotted as a function of j for various time differences r is shown in 



Fig. 15 for a soft chain and in Fig. 16 for a harmonic chain. For a given coupling constant k 
and delay time r, the correlation function peaks at the site i+j to which most of the energy 
originally at i has migrated. The change of the peak position with k indicates the velocity 
of the migration, and the height and width of the pulse reflect the dispersive dynamics. 
The following results are evident: 
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Figure 15: Dynamical energy correlation function C{j,T) for soft chains with 7 = 0.05 and 
ksT = 1.0. 
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Figure 16: Dynamical energy correlation function C(j, r) for harmonic chains with 7 = 0.05 
and ksT = 1.0. 
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• Increasing k in either soft or harmonic chains increases the velocity at which a fluc- 
tuation propagates. 

• The velocity for a given set of parameters is greater in the soft chain. 

• Dispersion is slower in the soft chain. 

However, as noted before, the differences between soft and harmonic chains at large k are 
fairly marginal. More dramatic differences in mobility features occur with anharmonic 
intermolecular potentials, a situation that will be presented elsewhere.@ 



6 Conclusions 

We have presented a fairly complete characterization of the thermal equilibrium behavior of 
oscillator chains with "diagonal anharmonicity," that is, chains with nonlinear on-site poten- 
tials and harmonic intersite potentials. Our particular interest lies in the characterization 
of possible spatial energy localization in such systems, and of the temporal persistence of 
such localization. 

The instantaneous localization of energy of a system in thermal equilibrium is a man- 
ifestation of the thermal fluctuations: it is an equilibrium property unrelated to system 
dynamics. We argued that not only do soft anharmonic chains have a higher total energy 
at a given temperature than do harmonic or hard chains, but also that thermal fluctuations 
are more pronounced in the soft anharmonic chains. This is a consequence of the fact that 
free energy maximization favors the occupation of phase space regions with a high density 
of states. The density of states increases with energy in a softening potential, so it is entrop- 
ically favorable for a few soft oscillators to have rather high energies. This in turn leads to 
greater spatial energy variability than in harmonic or hard chains, that is, soft chains have 
"hotter spots." The effect becomes more pronounced with increasing temperature. This 
entropic energy localization mechanism in soft chains is degraded as the harmonic intersite 
potential increases because the harmonic contributions become dominant over the local soft 
anharmonicity effects. 

In addition to the capacity for instantaneous localization of energy (which is greatest 
in soft chains), one is interested in the temporal degradation of a high energy fluctuation. 
That is, given a "local hot spot" (which is easier to find in soft chains, but nevertheless 
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does occur in harmonic and hard chains due to thermal fluctuations), how does such a 
fluctuation evolve in time? Such a fluctuation never grows spontaneously, nor does it 
persist indefinitely. Rather, it eventually degrades, either through dissipation into the bath 
or through dispersion along the chain. 

The rate of dissipation into the bath depends on the value of the dissipation parameter 
and also on the kinetic energy of the oscillators. If the dissipation parameter is small, 
this channel is of course slow for any chain. However, even if the dissipation parameter is 
large, dissipation can still be slow if the energy is not primarily in kinetic form. This is the 
case for soft chains provided the interatomic potential is weak (since otherwise the chain is 
essentially harmonic). In soft chains the energy is in potential form for a longer fraction of 
the time than in the other chains. As temperature is increased, this effect becomes more 
pronounced because ever softer portions of the potential become accessible, and the energy 
is stored as potential energy a greater fraction of the time. 

Thus an increase in temperature in weakly coupled soft chains leads not only to greater 
energy fluctuations but also to a slower decay of these fl,uctuations. 

Energy dispersion along the chain depends on the magnitude of the coupling constant 
and also on the relative oscillator displacements. If the coupling constant is small, this 
channel is slow for any chain. If it is large, then this channel can still be slow if relative 
displacements of neighboring oscillators are small. This is the case for the hard chain, 
where displacements are relatively small and don't change much with increasing energy. 
Furthermore, because in a hard oscillator the frequency increases with increasing energy, 
there is a frequency mismatch between a "hot" oscillator and its "colder" neighbors that 
further impedes energy transfer. This leads to greater persistence of local high-energy 
fluctuations with increasing temperature. 

Thus an increase in temperature in weakly dissipative hard chains leads not only to 
greater energy fluctuations but also to a slower decay of these fluctuations. 

The soft chain, on the other hand, increasingly loses its soft character as the interoscilla- 
tor energy transfer channel strengthens, and therefore both the landscape and the dynamical 
effects of anharmonicity quickly disappear as this coupling constant is increased. 

Finally, we showed that in harmonically coupled nonlinear chains (that is, in chains 
with "diagonal anharmonicity") in thermal equilibrium, high-energy fluctuation mobility 
does not occur beyond that which is observed in a harmonic chain. The situation might 
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be quite different if there is "nondiagonal anharmonicity" , that is, if the interoscillator 
interactions are anharmonic. Our results on these systems wih be presented elsewhere.0 

Further presentations will also deal with bistable "impurities" connected to chains of 
the types that we have considered hereS 
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